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Abstract. We present large-scale simulations of radial Eden clusters in three- 
dimensions and show that the growth exponent is in agreement with the value 
(3 = 0.242 accepted for the Kardar-Parisi-Zhang (KPZ) universality class. Our 
results refute a recent assertion proposing that radial Eden growth in d = 3 
belongs to a universality class distinct from KPZ. We associate the previously 
reported discrepancy to a slow convergence to the asymptotic limit. We also 
present the skewness and kurtosis in the roughening regime for flat geometry in 
2+1 dimensions. 



The Kardar-Parisi-Zhang (KPZ) universality class was introduced by the 
equation [T] 

f = »V 2 h + ±(Vhf + r,, (1) 

that describes the non-conservative evolution of an interface h(x,t) subjected to a 
white noise ?y. Many analytical [5J |3], numerical [H |H1 [7] and experimental [5] 
evidences show that radial (curved) and flat geometries have the same scaling 
exponents in 1+1 dimensions. 

The KPZ class is featured by an interface width w, defined as the standard 
deviation of the distance between the interface points and the initial seed, scaling 
in time as w ~ f 3 where is the growth exponent. In d = 1 + 1 dimensions 
the KPZ growth exponent is exactly known as f3 = 1/3 [1 . Eden model [5] is a 
standard KPZ system that generates radial clusters with a rough surface. Off-lattice 
simulations of Eden clusters immersed in d = 2 dimensions confirmed the agreement 
with the KPZ growth exponent |4j [7] . Recently, Kuennen and Wang [10] reported 
off-lattice simulations of radial Eden clusters in d — 3 dimensions, in which a growth 
exponent (3 w 0.1, much smaller than the accepted KPZ value for flat substrates in 
2+1 dimensions, (3kpz ~ 0.24 [TT], was observed. 

In order to verify the validity of this unexpected result, we have performed large- 
scale simulations of an off-lattice Eden model. The model is defined as follows. At 
t = 0, a single particle stuck to the origin is introduced. New particles are added to the 
cluster, one at a time, at random positions such that the new particle is adjacent to at 
least one particle in the cluster, but does not overlap any other added particle. The 
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Figure 1. (a) Surface of a radial Eden D cluster with 9 X 10 5 particles, (b) 
Surface obtained with on-lattice Eden D model after the growth of 2000 layers 
using a flat substrate of size 500 X 500. 

dynamics runs as follows. A particle of the cluster is randomly chosen and a random 
direction is uniformly chosen on a unitary sphere. A new particle is attached along 
this position if no overlap happens. Otherwise, the simulation proceeds to the next 
step with the choice of a new particle. Independently of the growth success, the time 
is incremented by At = 1/N p , where N p is the total number of particles in the cluster. 
This version was studied in two dimensions in Ref. [7] and we call it of Eden D, since 
there exist three other versions commonly called A, B and C [12]. Notice that we 
did not care about optimizations 1 13] to avoid uncontrolled sub-leading corrections. 
A typical surface of a small cluster is shown in figure [lja) . 

Figure [2] shows the interface width against time, averaged over 230 clusters. The 
largest sizes contain more than 10 s particles that correspond to a size much larger 
than 5 x 10 5 particles simulated in Ref. [TU]. Each sample takes about 8 hours of 
simulation in a Intel Xeon 2.53 GHz CPU and 7 GB of RAM memory. The effective 
exponent /3 e // obtained from the logarithmic derivative of w vs. t is shown in the 
inset of figure [2j The KPZ growth exponent is approached only in the limit of very 
large clusters. 

The model implementation in Ref. [TU] is slightly different from the Eden D used 
in the present work. There, when a randomly chosen particle has enough nearby 
space to add a new particle, the growth always happens. This version is know in 
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Figure 2. Main Plot: Interface width against time for off-lattice simulation of 
radial Eden D in d = 3 dimensions. Dashed line has a slope 0.24 as a guide to the 
eyes. Inset: Local slope in plot In id vs. Int. Dashed line is the KPZ exponent 
P = 0.24. 



the literature as Eden B [T3]. In Eden D, frustrated attempts are allowed, which, 
consequently, amplifies the noise and generates larger interface fluctuations when 
compared to Eden B. 

In order to shed light on this issue, we have also performed flat simulations of both 
Eden models on square lattices with periodic boundary conditions along directions x 
and y. The initial condition consists in all sites occupied for z < 0. In Eden B, 
particles in the interface (an occupied site with at least one empty nearest neighbour 
(NN)) are randomly chosen and one of their empty NNs, also chosen at random, is 
occupied. Eden D is equal to Eden B except that any NN (occupied or not) can be 
selected. If an occupied neighbour is chosen, the simulation runs to the next step. In 
both models the time step is At — 1/N S , where N s is the number of interface sites. A 
typical surface for Eden D model is shown in figure [ijb) . 

Figure [3] shows the interface width against time for both flat models. Averages 
were performed over 400 independent samples. The inset shows the corresponding 
effective growth exponent. While the growth exponent in Eden D simulations 
converges quickly to a value close to the KPZ universality class, in Eden B, an 
exponent increasing from /3 > 0.1 but not exhibiting a stationary power regime is 
observed. These results show that Eden B model has very strong corrections to the 
scaling and that the KPZ exponent must be observed only in exceedingly long time 
simulations. It is worth noticing that the interface width in Eden B is very small 
(w < 2 lattice unities in figure [3] and w < 1.5 particle diameters in Ref. QZ3), which 
hinders the observation of the asymptotic regime. 

Recent advances in the understanding of the KPZ universality class in 1+1 
dimensions have been based on universal distributions of the interface fluctuations 
during the roughening (transient) regime [2]. A basic characterization of the 



Alves and Ferreira, Eden clusters in three- dimensions and the KPZ class 
.1 



10 



10 



0.5 
0.4 

to 0.3 

o 

02- 0.2 
0.1 





10 







T 



T 



T 



°o°°ocW 



10 



_l , L 

4 6 

lnt 



10 
t 



EdenB 
EdenD 



10 



10 



10 



Figure 3. Main Plot: Interface width against time for lattice simulations of 
Eden B and D on flat substrates of size 1000 X 1000. The dashed line has a slope 
0.24 as a guide to the eyes. Inset: Local slope in plots In mi vs. Int. Dashed lines 
represent the KPZ exponent /3 = 0.24 and the estimate of Ref. 1101 /3 = 0.1. 



distributions can be obtained 
dimensionless cumulant ratios 

(h 3 )c 



via skewness S and kurtosis K defined as the 



and 



S = 



K = 



-3, 



(h 2 )l 

where (X n ) c denotes the nth cumulant of X. In our simulations, these quantities did 
not exhibit stationary values neither for radial Eden D (Fig. [4]) nor for Eden B on 
flat substrates. For Eden D on flat substrates, the skewness and kurtosis converge to 
S = 0.411(2) and K — 0.34(1), respectively. Figure [4] shows the skewness and kurtosis 
against time for radial and flat Eden D simulations. It is important mentioning 
that cumulant ratios for Eden D in flat simulations have finite size effects. So, it 
was necessary to run shorter times and systems larger than those used in figure [3] 
to obtain an accurate estimate of the cumulant ratios. The skewness and kurtosis 
obtained Eden D agree pretty well with the ratios obtained for other KPZ models in 
2+1 dimensions [T5]- A detailed analysis of the height distributions for several KPZ 
models in three-dimensional lattices will be reported elsewhere |15) . 

In conclusion, we have shown that the radial Eden growth in three-dimensions 
has a growth exponent in agreement with the KPZ universality class, at odds with a 
previous report [10] but in consonance with the common belief. Our results suggest 
that the low exponent observed in Ref. |10) comes from the small simulated sizes and 
strong corrections to the scaling. In addition, we have presented the stationary values 
of the skewness and kurtosis during the roughening regime for flat simulations. 
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Figure 4. Main plot: Skewness against time for radial (off-lattice) and flat 
(on-latticc) simulations of the Eden D model. Flat simulations were performed 
in lattices with 2200 X 2200 sites. Inset: Corresponding kurtosis against time. 
Dashed lines represent the estimated cumulant ratios. 
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